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Abstract 

The quantum kinetic equation used in the study of weak turbulence 
is reconsidered in the context of a theory with a generic quartic inter- 
action. The expectation value of the time derivative of the mode num- 
ber operators is computed in a perturbation expansion which places 
the large diagonal component of the quartic term in the unperturbed 
Hamiltonian. Although one is not perturbing around a free field the- 
ory the calculation is easily tractable owing to the fact that the un- 
perturbed Hamiltonian can be written solely in terms of the mode 
number operators. 
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1 Introduction 



In one approach to the statistical description of weak turbulence, a central 
role is played by the kinetic wave equation |], ||. This equation for the time 
derivative of the mode numbers has been derived for both classical and quan- 
tum systems in a perturbation series by expanding about a free field (har- 
monic oscillator) theory. In this paper, we will reconsider this derivation for a 
quantum mechanical system whose Hamiltonian is a sum of generic quadratic 
and quartic terms. Our perturbation expansion will perturb around an op- 
erator which contains the diagonal component of the quartic term together 
with the usual quadratic, or free field, component. Since one is interested 
in the expectation values of fields between states with large mode numbers, 
it is sensible to include as much of these in the unperturbed Hamiltonian 
as the calculation permits. We do not need to assume that the coupling to 
all the quartic terms in the Hamiltonian is small; the diagonal part can be 
arbitrarily large in our approach. 

We begin in the following section with a precise statement of the theory 
under consideration and an encapsulation of the interaction picture used to 
carry out the derivation of the quantum kinetic equation. At this order, the 
largest terms are cubic in the mode numbers and these we calculate explic- 
itly. We will find some additional terms not present in another derivation 
of the quantum kinetic equation [p]], and then go on to consider under what 
conditions one might expect those corrections to be small. A discussion of 
stationary solutions to the kinetic equation then follows. 



2 The Quantum Kinetic Equation 

Let us consider a quantum mechanical system based on the Hamiltonian, 

H = ^2 u k a{ a k + ^2 T klk2kaki a kl a\ 2 a ks a k4 , (1) 

k ki,—,kt 

which contains generic quadratic and quartic terms; the number d of spatial 
dimensions in which the system evolves is arbitrary. The free field oscilla- 
tor energies uj k are assumed given and constitute part of the specification of 
the system. The function T klk2k3k4 includes, by definition, the momentum 
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conserving factor 5 kl+k , 2tk . i+ki (both the vector ki and the 5-function should 
be understood as d- dimensional quantities). Beyond the implicit symmetry 
properties which are that T klk2kzkA is symmetric in the first two and last two 
indices, and that under complex conjugation T klk2k3k4 = T k3k4klk2 , this coef- 
ficient may contain further momentum dependence which we will otherwise 
not restrict in the derivation of the kinetic equation. As usual, a k and a k in 
eqn. ([!]) denote Bose creation and annihilation operators, and they obey the 
commutation relation, 

[a k ,a]} = S k ,i ■ (2) 



It is convenient to use the number operator n k = a k a k for each mode in 
our system. The states which diagonalize these number operators satisfy [§]: 
h k \n k ) = n k \n k ), a k \n k ) = ^/n^\n k -l), and a{ \n k ) = y/n k + 1 \n k + l), and 
are clearly labelled by their eigenvalues. 

A perturbation series (see || for a thorough exposition) begins by split- 
ting the Hamiltonian into an unperturbed component H and a "small" com- 
ponent Hi. In our approach, we will place the diagonal of the quartic com- 
ponent in the unperturbed sector, so that H = H + Hi with, 

H = ( u k ~ 2 T k) hk + Z T ki ™* hi (3) 

k k,l 

Hi = ^ ^fcifc 2 fc 3 fc 4 a fci a fc 2 flfc 3 a kA 5 

and where we have introduced the notation T k = T kkkk , T k i = T k i k i, and 

rpi ( T klk2k3k4 if ki ^ k 3 or k 4 , . 

k!k 2 k 3 k 4 | o otherwise ' 1 ' 

We can express this equivalently as, 

Tk^kski = ( 1 _ <W 3 <W 4 ~ fikiki 4 2 fc 3 + <W 2 faik 3 <W 4 ) T klk2k3k4 . (5) 

It is important to emphasize that we need not assume that the coefficient 
functions Tjy be small, as they are part of the unperturbed Hamiltonian. 
The validity of the perturbation expansion depends, however, that T kik2kak4 
be "small" relative to the mode numbers. 
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In the Heisenberg representation of quantum mechanics, operators are 
time dependent while the states are time independent. Given some operator 
A, its time evolution as a Heisenberg operator A H (t) satisfies 

jA H (t)=i[H,A H (t)] , (6) 

and one may equivalently write this as A B {t) = exp[iHt] Ah{0) exp[— iHt]. 
Any expectation value in this representation therefore satisfies, 

j t (*i|A H (t)|* 2 > = (*x| j t A H {t) |* 2 > = {^[H.A^t)] |* 2 > , (7) 

since the states |\l/ a ) are independent of t. We are interested in the case 
A = h k for large t, and we will compute, 

lim (*\i[H,n kB (t)] |*> , (8) 

t — >oo 

for some state \^) which we will later specify. This is the precise meaning 
we ascribe to the time derivative of the mode number appearing in other 
presentations |], of the kinetic wave equation. 

For the purposes of perturbation theory, one moves to the interaction 
picture where the following relations hold, 

(*x| H (t) |* 2 ) = Oj(t) |* 2 (t)> (9) 

|* a (t)) = exp[iH t) exp[— z H (t — t')] exp[-iH t'] |* a (t')) 

0/(t) = exp[t H (t - t')] 0/(0 exp[-i F (t - 0] • 

for all operators O and all states \^ a )- Interestingly, the time evolution of 
the operators a kI and a^i in this model is very simple in spite of the fact 
that they do not evolve via a free field Hamiltonian. It is not difficult to first 
show that 

[ojb, H ] = (uj k + T k i hi) a k , (10) 
l 

and using this one can quickly prove, 

a>ki(t) = exp[iH t] a kI (0) exp[—iH Q t] (11) 
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= exp[-2 1 {u k + 4 T kt hi)} a kI (0) 

i 

= a kI (0) exp[-it(uj k -AT k + Aj2T H ni)] . 

i 

The combination of operators in the second line of @ is conveniently 
denoted by, 

U(t,t') = exp[iH t] exp[-i H (t - t')} exp[-iH Q t') (12) 

= 1 -i C dTH u (r) U(r,t') . 
Jt' 

To lowest order in the interaction Hi, we just set U(T,t r ) = 1 on the right 
hand side of eqn. flT2|). Our goal is to compute, 

(j t n k ) = hm U\t,-t) Oj(t) U(t,-t) \*(-t)) , (13) 

where O = i[Hi,h k ] . We should emphasize that the states on both sides of 
this expectation value are m-states in the sense of scattering theory as both 
are at — oo. 

Noticing that H commutes with h k , the computation of eqn. ( |i~3"D to the 
lowest order in perturbation theory reduces to, 

lim (tf(-t)| i[H u (t),h k } + f [[H u {t\n k lH u {r)\dr \*{-t)) . (14) 

We assume that the state |^ r (— oo)) is an eigenstate of the number operators 
n k . It is not difficult to show that the first order term in flT4]) does not 
contribute; to see this one simply computes, 

[Hi, h k }=2 (T'k^ksk a ki a k 2 Gfc 3 a k - T kl kk 2 k 3 a ki a k a k 2 a k 3 ) ■ (15) 

kik 2 k 3 

When we take this between the same states, there is no way to pair the 
creation and annihilation operators as T' is off-diagonal, and hence the first 
term in (|14|) manifestly vanishes. However, we should point out that this 
term would also vanish even if we had not subtracted out the diagonal, and 
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had instead made the normal perturbative expansion around the quadratic 
term in H; the cancellation would then involve a mixing of both terms in 
eqn. flI5p. 



One of the ingredients needed in this computation is the expectation 
value, 

(4i 4 2 a k 3 a k4 «?! a h a h a h) > ( 16 ) 

where k\ or k 2 7^ k 3 or fc 4 , and l\ or l 2 7^ ^3 or Z 4 . The states on both 
sides are identical and are eigenstates of all the number operators. A simple 
computation yields, 

(n h + 1) n h [ (8 hk3 8 hk/k + 8 hki 8 hk3 ) (n h + 1) - 8 hh 8 hk3 8 hki n h ] • 

[ {8i 3 k! 8i 4 k 2 + 8i 3 k 2 <W) n>i 4 - 8i 3 i 4 5i 3kl 8i ik2 (n k + 1) ] . (17) 

We should note that it is not sufficient for the purposes of our calculation, 
even working in the large ni limit, to keep only the most dominant terms in 
the above expression which are fourth order in these mode numbers. When 
this expression is used in our calculation, we will see that the fourth order 
terms cancel and the next to leading terms remain. For this reason we have 
been careful to take account of the possibility that li equals l 2 , and so on, in 
this expectation value; no assumptions (e.g. random phase approximation; 
see [Q]) have been made in obtaining the expression in flTT|). 

It is straightforward, though rather tedious, to assemble the above pieces, 
and the tree level expression for the quantity in ( |13D is found to be, 



(j^k) = 87r J2 \ T Lk 2 k 3 k\ 2 8{kx,k 2 ,k 3 ,k) [ s 3 (k u k 2 , k 3 , k) (18) 

ki,k 2 ,k3 

+s 2 (h, k 2 , h, k) + sx(ki, k 2 , k 3 , k) ] , 

with the functions s a (ki, k 2 , k 3 , k) given by, 

S3 = 4 ( n kl n k2 n k3 + n kl n k2 n k - n kl n k3 n k - n k2 n kj n k ) (19) 
- 2 8 klk2 ( n kl n k2 n k3 + n kl n k2 n k ) + 2 8 kak ( n kl n H n k + n k2 n ks n k ) 

S2 = 4 ( n kl n k2 - n ks n k ) - 2 5 klk2 ( n kl n k2 + n kl n k + n kl n k , A ) 
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+ 2 5 k3k ( n kl n k3 + n k , 3 n k + n k2 n ka ) 

Si = —2 5 kl k 2 n kl + 2 5 k , ik n ka ■ 

For large values of the mode numbers, the S3 term which is cubic in those 
variables will dominate. The energy conserving delta function S(ki, fc 2 , k 3 , k), 
which arose from an integration over r, is given by, 

5(u kl + u k2 - u k:j - ui k + 4 (Tkii + T k2 i - T k3 i - T kl ) n L (20) 

1 

+ 4 (Tfe 3 + T k + T klk2 + T kak — T klk — T k2k — T klk , 3 — T k2ka ) ) . 

It is obtained by moving the r operator dependence entirely to either the 
left or the right where it becomes a normal function of r after acting on 
the states, and then using the representation, f_ ™ e tTX dr = 2n S(x). This 
rather asymmetric looking expression can be recast as, 

S( {uo kl - 2T kl ) + {u k2 - 2T k2 ) - {u k3 - 2T k3 ) - (u k - 2T k ) (21) 

+ 2^2 {T kl i + T k2 i - T k3 i - T k {)(2ni + 5 kll + 5 k2 i - 5 k3 i - 5 k i) ) • 
1 

If we define an effective energy per state by 

£ m = u m - 2T m + 2 T m i ( 2 rii + 5 kll + 5 k2l - 5 k . 3 i - 5 kl ) , (22) 

1 

then the delta function can be written most simply as, 

5( e kl + e k2 - e ka - e k ) . (23) 

The exact expression for e m given above is greatly simplified in the large rii 
situation in which we work, so that 

e m ss w m + 4^ T ml n { . (24) 
1 

Let us emphasize that the leading order terms which are cubic in the 
mode numbers in flTB|) do not fully agree with the expressions found in [fll 0]; 



6 



our leading order terms can be written as, 



(-r.hk) = 8tt \ T U 2 k 3 k\ 2 S(e kl +£k 2 ~e k3 - e k ) • (25) 

ki,k 2 ,k 3 

r a, 1 1 1 1 ^ 

n kl n k2 n k3 n k { 4 ( 1 ) 

n k3 n k n kl n k , 2 

-2 8 klk2 ( — + —) + 2 5 k3k ( — + —)} , 
n kz n k n kl n k2 

with e m given by eqn. (0). One clear difference in our result concerns the 
last two terms in (|25|) which are not present in p] ; these would still be there 
in our analysis even if we had perturbed around the quadratic term in H . In 
fact, further diagonal terms would presumeably be present as well, but these 
we have accounted for by incorporating them into H . This discrepancy will 
be discussed in the following section. The other clear difference is due to our 
different perturbation expansion which results in the 5-function involving the 
effective energy e k , rather than uj k . 



3 Stationary Solutions 

The analysis of stationary solutions to the kinetic equation in the usual per- 
turbative version can be found in [p], |2| . In that analysis, the extra terms we 
found in (|25|) are not considered, and they analyze the condition, 

= l^fcife 2 fc 3 fc| 2 8(io kl + u k2 — uj ka — to k ) ■ (26) 

ki, k 2,k2, 

r 1 1 1 1 i 

n kl n k2 n k3 n k { 1 } . 

n k:j n k n kl n k2 

The question therefore arises as to whether the other terms in eqn. (|25|), 
which are also cubic in the mode numbers, are small relative to the others. 
In this regard, it is noteworthy that the additional terms involve one less sum 
over momentum space. So if the terms in the summation are large over some 
reasonable domain in momentum space, then those factors will be suppressed 
by roughly the volume of that domain. Our calculation thus far has assumed 
that at least some of the mode numbers are large compared to unity and 
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neglecting these additional two terms amounts to some kind of additional 
condition such as the one just suggested. One might instead entertain a 
more restricted class of T' kik2hiki coefficient such that the additional terms 
vanish identically; perhaps this might be natural in the context of vortex 
dynamics. For example, it could contain the factor \k\ — k^\° ■ \k 3 — k^, 
for a > 0. We will not consider this issue further here, and will proceed 
to consider the stationary solutions to eqn. fl25|) assuming that the last two 
terms can be neglected. 

Following the analysis in the usual perturbative expansion [jl], j2|, one 
solution for the occupation numbers n m is given by, 

T , , 



where T and \i are constants; this is the large T limit of the usual thermody- 
namic distribution (exp[(ek + fi)/T] — || for noninteracting bosons. The 
difference in our case is simply that the effective energy e m enters rather than 
uj m . Given the precise form of this effective energy, we see that eqn. ([27]) is 
in fact a self-contained integral equation for n m , but we will not analyze it 
further here. 

While the preceeding analysis can equally well be carried out in terms of 
discrete momenta and sums (which we have done), and continuous variables 
with integrals, the examination of the Kolmogorov solutions requires the later 
setting. The only care in going over to integral expressions is in correctly 
treating 5-function factors. It is generally the case that the perturbative 
expansion we have considered will contain a factor of 5(0) which must be 
factored out and discarded, and it is therefore convenient now to take the 
ik 2 k 3 k 4 coefficient without the momentum conserving ^-function, 5^{ki + 
^2 ~ k 3 — k), we previously included so that the integral equivalent of eqn. 
( |25D only has a single momentum conserving ^-function. 

One way to establish the Kolmogorov solutions is to follow the presen- 
tation in [|TJ substituting the effective energy e m for the free field energy 
uj m . This entails a number of assumptions. We will assume that the theory 
has rotational symmetry which implies, in particular, that the mode number 
rik = n(k) and the effective energy e fc = e(k) only depend on the magnitude of 
the vector k. We further suppose that the scaling properties, e(X k) = X a e(k) 
and T'(A k\, X k 2 , X k 3 , X k±) = X 13 T'(ki, k 2 , k%, k^) are satisfied, and that the 
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functional relation e(k) is invertible. Without loss of generality, we can take 
e(0) = by adjusting, say, the u>(0) coefficient if necessary. Given these as- 
sumptions, we will show that a solution exists of the form n(k) = n{e) = e~ x 
(repeated use of the symbol n for two different functions should cause no con- 
fusion; we will always regard the mode numbers as functions of the effective 
energy in the following). 



Let us begin by integrating over angles in (P5|); the volume element is 
d d ki = kf~ x dki d£li and it should be clear from context whether we use the 
symbol kj to denote a vector or its magnitude. Using the assumption that 
we can invert the relationship e(k), we can change variables in the remaining 
integrals from k{ to £j. It is convenient to first define, 

ms u e»e.,e) = l£^||-' (28, 



/ K k2 ksk\ 2 (h + k 2 — k 3 — k) dttt dtt 2 dQ 3 



where we are integrating over three sets of angular variables. Let us note 
that U shares the same symmetry properties as the coefficient T' under per- 
mutations of its arguments. Our task now is to find solutions to, 







poo 

l de 1 de 2 de 3 U(e 1 ,e 2 ,e 3 ,e) 6(e 1 + e 2 - e 3 - e) (29) 
Jo 

r 1 1 1 1 -i 

n kl n k2 n k3 n k { 1 } . 

n ka n k n kl n k2 

It is straightforward to work out the scaling properties of U which are 
implied by the assumptions. By simply scaling each of the momenta on both 
sides of the defining equation, one quickly finds the relation, 

L/"(Aei,Ae a ,Ae3,Ae) = A 7 U(e x , e 2 , e 3 , e) (30) 

3d + 2(3 A 
7 = - 4 . 

a 



Returning to the analysis of eqn. , one easily uses the 5-function to 
carry out the e 3 integral; what remains is a region D of the (ei, e 2 )-plane. 
This region is not the entire first quadrant since we must satisfy the condition, 

e 3 = £i +e 2 -e > . (31) 
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Divide D into four sectors as follows, 

Da = {(e 1 ,e 2 )ED | e x < e, e 2 < e} (32) 

D 2 = {(£ 1 ,e 2 )eD | e 1 >e, e 2 > e} 

D 3 = {(£ 1) £ 2 )eD | ei < e, £ 2 >e} 

L> 4 = {(ei,e 2 )eD | £i>£, £ 2 < , 

and perform the Zakharov transformations []2], f|] to map D 2 , D3, and D4 
onto D\. Those transformations respectively take the form, 

D2 : e t = , f] e 2 = , ff , (33) 

e\ + e 2 — e e\ + e 2 — e 

e(e[ + e' 2 -e) e 2 
■ £1 — -, £2 — —r , 

_ £ 2 e(ei + 4 - e) 

JJ4 . £1 — — e 2 







Using these transformations, and the ansatz n(e) = C e x , one finds that 
eqn. fl29D becomes, 

/ deide 2 U(e 1 ,e 2 ,e l + e 2 -e,e) [e 1 e 2 (e 1 + e 2 -e)e}' x (34) 
{ {si + e 2 - ef + e x - e\ - e x 2 } 

{i+( £i+ ; 2 ~7 -(f) y -(^} , 

with y = 3x — 7 — 3. We can satisfy this relation at the four points x = 0, x = 
1, y = 0, or y = 1. The first two cases are simply limits of the thermodynamic 
distribution previously considered, while the last two possibilities correspond 
to the Kolmogorov solutions, 



n 



(i) 



( £ )= Cl e-^ 3 , n^(e) = C 2 s-^ 3 . (35) 



Although the form of these solutions is the same as in the usual perturbative 
analysis, they differ fundamentally in that they are in fact integral equations 
for n(s), owing to eqn. we will not pursue an analysis of them here. 
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4 Concluding Remarks 



In this paper we have established a different perturbative expansion which 
leads to a quantum kinetic equation, similar to the usual one obtained by 
expanding about a quadratic Hamiltonian. By placing the diagonal quartic 
terms into the unperturbed Hamiltonian, we have generalized that equation 
and identified the effective state energy which essentially takes the place of 
the free field energy in those earlier treatments. The framework we establish 
applies as well to the usual perturbative situation and can be considered as 
an alternative to other derivations. We have also highlighted some subtleties 
that arise there in the leading order terms cubic in the mode numbers and 
have offered an argument why one might expect those additional corrections 
to be small. 
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